Systematic analysis of RNA-binding proteins identifies targetable therapeutic vulnerabilities in osteosarcoma.

Osteosarcoma is the most common primary malignant bone tumor with a strong tendency to metastasize, limiting the prognosis of affected patients. Genomic, epigenomic and transcriptomic analyses have demonstrated the exquisite molecular complexity of this tumor, but have not sufficiently defined the underlying mechanisms or identified promising therapeutic targets. To systematically explore RNA-protein interactions relevant to OS, we define the RNA interactomes together with the full proteome and the transcriptome of cells from five malignant bone tumors (four osteosarcomata and one malignant giant cell tumor of the bone) and from normal mesenchymal stem cells and osteoblasts. These analyses uncover both systematic changes of the RNA-binding activities of defined RNA-binding proteins common to all osteosarcomata and individual alterations that are observed in only a subset of tumors. Functional analyses reveal a particular vulnerability of these tumors to translation inhibition and a positive feedback loop involving the RBP IGF2BP3 and the transcription factor Myc which affects cellular translation and OS cell viability. Our results thus provide insight into potentially clinically relevant RNA-binding protein-dependent mechanisms of osteosarcoma.

Evolutionary mode and timing of dissemination of high-grade serous carcinomas.

Dissemination within the peritoneal cavity is a main determinant of poor patient outcomes from high-grade serous carcinomas (HGSCs). The dissemination process is poorly understood from a cancer evolutionary perspective. We reconstructed the evolutionary trajectories across a median of 5 tumor sites and regions from each of 23 patients based on deep whole-exome sequencing. Polyclonal cancer origin was detected in 1 patient. Ovarian tumors had more complex subclonal architectures than other intraperitoneal tumors in each patient, which indicated that tumors developed earlier in the ovaries. Three common modes of dissemination were identified, including monoclonal or polyclonal dissemination of monophyletic (linear) or polyphyletic (branched) subclones. Mutation profiles of initial or disseminated clones varied greatly among cancers, but recurrent mutations were found in 7 cancer-critical genes, including TP53, BRCA1, BRCA2, and DNMT3A, and in the PI3K/AKT1 pathway. Disseminated clones developed late in the evolutionary trajectory models of most cancers, in particular in cancers with DNA damage repair deficiency. Polyclonal dissemination was predicted to occur predominantly as a single and rapid wave, but chemotherapy exposure was associated with higher genomic diversity of disseminated clones. In conclusion, we described three common evolutionary dissemination modes across HGSCs and proposed factors associated with dissemination diversity.

The interplay between neoantigens and immune cells in sarcomas treated with checkpoint inhibition.

Introduction: Sarcomas are comprised of diverse bone and connective tissue tumors with few effective therapeutic options for locally advanced unresectable and/or metastatic disease. Recent advances in immunotherapy, in particular immune checkpoint inhibition (ICI), have shown promising outcomes in several cancer indications. Unfortunately, ICI therapy has provided only modest clinical responses and seems moderately effective in a subset of the diverse subtypes. Methods: To explore the immune parameters governing ICI therapy resistance or immune escape, we performed whole exome sequencing (WES) on tumors and their matched normal blood, in addition to RNA-seq from tumors of 31 sarcoma patients treated with pembrolizumab. We used advanced computational methods to investigate key immune properties, such as neoantigens and immune cell composition in the tumor microenvironment (TME). Results: A multifactorial analysis suggested that expression of high quality neoantigens in the context of specific immune cells in the TME are key prognostic markers of progression-free survival (PFS). The presence of several types of immune cells, including T cells, B cells and macrophages, in the TME were associated with improved PFS. Importantly, we also found the presence of both CD8+ T cells and neoantigens together was associated with improved survival compared to the presence of CD8+ T cells or neoantigens alone. Interestingly, this trend was not identified with the combined presence of CD8+ T cells and TMB; suggesting that a combined CD8+ T cell and neoantigen effect on PFS was important. Discussion: The outcome of this study may inform future trials that may lead to improved outcomes for sarcoma patients treated with ICI.

Chromosomal instability and a deregulated cell cycle are intrinsic features of high-risk gastrointestinal stromal tumours with a metastatic potential.

Patients with localised, high-risk gastrointestinal stromal tumours (GIST) benefit from adjuvant imatinib treatment. Still, approximately 40% of patients relapse within 3 years after adjuvant therapy and the clinical and histopathological features currently used for risk classification cannot precisely predict poor outcomes after standard treatment. This study aimed to identify genomic and transcriptomic profiles that could be associated with disease relapse and thus a more aggressive phenotype. Using a multi-omics approach, we analysed a cohort of primary tumours from patients with untreated, resectable high-risk GISTs. We compared patients who developed metastatic disease within 3 years after finishing adjuvant imatinib treatment and patients without disease relapse after more than 5 years of follow-up. Combining genomics and transcriptomics data, we identified somatic mutations and deregulated mRNA and miRNA genes intrinsic to each group. Our study shows that increased chromosomal instability (CIN), including chromothripsis and deregulated kinetochore and cell cycle signalling, separates high-risk samples according to metastatic potential. The increased CIN seems to be an intrinsic feature for tumours that metastasise and should be further validated as a novel prognostic biomarker for high-risk GIST.

Heritable defects in telomere and mitotic function selectively predispose to sarcomas.

Cancer genetics has to date focused on epithelial malignancies, identifying multiple histotype-specific pathways underlying cancer susceptibility. Sarcomas are rare malignancies predominantly derived from embryonic mesoderm. To identify pathways specific to mesenchymal cancers, we performed whole-genome germline sequencing on 1644 sporadic cases and 3205 matched healthy elderly controls. Using an extreme phenotype design, a combined rare-variant burden and ontologic analysis identified two sarcoma-specific pathways involved in mitotic and telomere functions. Variants in centrosome genes are linked to malignant peripheral nerve sheath and gastrointestinal stromal tumors, whereas heritable defects in the shelterin complex link susceptibility to sarcoma, melanoma, and thyroid cancers. These studies indicate a specific role for heritable defects in mitotic and telomere biology in risk of sarcomas.

Mutational analysis and protein profiling predict drug sensitivity in multiple myeloma cell lines.

Introduction: Multiple myeloma (MM) is a heterogeneous disease where cancer-driver mutations and aberrant signaling may lead to disease progression and drug resistance. Drug responses vary greatly, and there is an unmet need for biomarkers that can guide precision cancer medicine in this disease. Methods: To identify potential predictors of drug sensitivity, we applied integrated data from drug sensitivity screening, mutational analysis and functional signaling pathway profiling in 9 cell line models of MM. We studied the sensitivity to 33 targeted drugs and their association with the mutational status of cancer-driver genes and activity level of signaling proteins. Results: We found that sensitivity to mitogen-activated protein kinase kinase 1 (MEK1) and phosphatidylinositol-3 kinase (PI3K) inhibitors correlated with mutations in NRAS/KRAS, and PI3K family genes, respectively. Phosphorylation status of MEK1 and protein kinase B (AKT) correlated with sensitivity to MEK and PI3K inhibition, respectively. In addition, we found that enhanced phosphorylation of proteins, including Tank-binding kinase 1 (TBK1), as well as high expression of B cell lymphoma 2 (Bcl-2), correlated with low sensitivity to MEK inhibitors. Discussion: Taken together, this study shows that mutational status and signaling protein profiling might be used in further studies to predict drug sensitivities and identify resistance markers in MM.

Strategies to inhibit FGFR4 V550L-driven rhabdomyosarcoma.

BACKGROUND: Rhabdomyosarcoma (RMS) is a paediatric cancer driven either by fusion proteins (e.g., PAX3-FOXO1) or by mutations in key signalling molecules (e.g., RAS or FGFR4). Despite the latter providing opportunities for precision medicine approaches in RMS, there are currently no such treatments implemented in the clinic. METHODS: We evaluated biologic properties and targeting strategies for the FGFR4 V550L activating mutation in RMS559 cells, which have a high allelic fraction of this mutation and are oncogenically dependent on FGFR4 signalling. Signalling and trafficking of FGFR4 V550L were characterised by confocal microscopy and proteomics. Drug effects were determined by live-cell imaging, MTS assay, and in a mouse model. RESULTS: Among recently developed FGFR4-specific inhibitors, FGF401 inhibited FGFR4 V550L-dependent signalling and cell proliferation at low nanomolar concentrations. Two other FGFR4 inhibitors, BLU9931 and H3B6527, lacked potent activity against FGFR4 V550L. Alternate targeting strategies were identified by RMS559 phosphoproteomic analyses, demonstrating that RAS/MAPK and PI3K/AKT are essential druggable pathways downstream of FGFR4 V550L. Furthermore, we found that FGFR4 V550L is HSP90-dependent, and HSP90 inhibitors efficiently impeded RMS559 proliferation. In a RMS559 mouse xenograft model, the pan-FGFR inhibitor, LY2874455, did not efficiently inhibit growth, whereas FGF401 potently abrogated growth. CONCLUSIONS: Our results pave the way for precision medicine approaches against FGFR4 V550L-driven RMS.

Validation of a novel risk score to predict early and late recurrence in solitary fibrous tumour.

BACKGROUND: Current risk models in solitary fibrous tumour (SFT) were developed using cohorts with short follow-up and cannot reliably identify low-risk patients. We recently developed a novel risk model (G-score) to account for both early and late recurrences. Here, we aimed to validate the G-score in a large international cohort with long-term follow-up. METHODS: Data were collected from nine sarcoma referral centres worldwide. Recurrence-free interval (RFi) was the primary endpoint. RESULTS: The cohort comprised 318 patients with localised extrameningeal SFTs. Disease recurrence occurred in 96 patients (33%). The estimated 5-year RFi rate was 72%, and the 10-year RFi rate was 52%. G-score precisely predicted recurrence risk with estimated 10-year RFi rate of 84% in low risk, 54% in intermediate risk and 36% in high risk (p < 0.001; C-index 0.691). The mDemicco (p < 0.001; C-index 0.749) and SalasOS (p < 0.001; C-index 0.674) models also predicted RFi but identified low-risk patients less accurate with 10-year RFi rates of 72% and 70%, respectively. CONCLUSIONS: G-score is a highly significant predictor of early and late recurrence in SFT and is superior to other models to predict patients at low risk of relapse. A less intensive follow-up schedule could be considered for patients at low recurrence risk according to G-score.

Clonal evolution in primary breast cancers under sequential epirubicin and docetaxel monotherapy.

BACKGROUND: Subclonal evolution during primary breast cancer treatment is largely unexplored. We aimed to assess the dynamic changes in subclonal composition of treatment-naïve breast cancers during neoadjuvant chemotherapy. METHODS: We performed whole exome sequencing of tumor biopsies collected before, at therapy switch, and after treatment with sequential epirubicin and docetaxel monotherapy in 51 out of 109 patients with primary breast cancer, who were included in a prospectively registered, neoadjuvant single-arm phase II trial. RESULTS: There was a profound and differential redistribution of subclones during epirubicin and docetaxel treatment, regardless of therapy response. While truncal mutations and main subclones persisted, smaller subclones frequently appeared or disappeared. Reassessment of raw data, beyond formal mutation calling, indicated that the majority of subclones seemingly appearing during treatment were in fact present in pretreatment breast cancers, below conventional detection limits. Likewise, subclones which seemingly disappeared were still present, below detection limits, in most cases where tumor tissue remained. Tumor mutational burden (TMB) dropped during neoadjuvant therapy, and copy number analysis demonstrated specific genomic regions to be systematically lost or gained for each of the two chemotherapeutics. CONCLUSIONS: Sequential epirubicin and docetaxel monotherapy caused profound redistribution of smaller subclones in primary breast cancer, while early truncal mutations and major subclones generally persisted through treatment. TRIAL REGISTRATION: ClinicalTrials.gov, NCT00496795 , registered on July 4, 2007.

Clonal evolution after treatment pressure in multiple myeloma: heterogenous genomic aberrations and transcriptomic convergence.

We investigated genomic and transcriptomic changes in paired tumor samples of 29 in-house multiple myeloma (MM) patients and 28 patients from the MMRF CoMMpass study before and after treatment. A change in clonal composition was found in 46/57 (82%) of patients, and single-nucleotide variants (SNVs) increased from median 67 to 86. The highest increase in prevalence of genetic aberrations was found in RAS genes (60% to 72%), amp1q21 (18% to 35%), and TP53 (9% to 18%). The SBS-MM1 mutation signature was detected both in patients receiving high and low dose melphalan. A total of 2589 genes were differentially expressed between early and late samples (FDR < 0.05). Gene set enrichment analysis (GSEA) showed increased expression of E2F, MYC, and glycolysis pathways and a decreased expression in TNF-NFkB and TGFbeta pathways in late compared to early stage. Single sample GSEA (ssGSEA) scores of differentially expressed pathways revealed that these changes were most evident in end-stage disease. Increased expression of several potentially targetable genes was found at late disease stages, including cancer-testis antigens, XPO1 and ABC transporters. Our study demonstrates a transcriptomic convergence of pathways supporting increased proliferation and metabolism during disease progression in MM.

miR-486-5p expression is regulated by DNA methylation in osteosarcoma.

BACKGROUND: Osteosarcoma is the most common primary malignant tumour of bone occurring in children and young adolescents and is characterised by complex genetic and epigenetic changes. The miRNA miR-486-5p has been shown to be downregulated in osteosarcoma and in cancer in general. RESULTS: To investigate if the mir-486 locus is epigenetically regulated, we integrated DNA methylation and miR-486-5p expression data using cohorts of osteosarcoma cell lines and patient samples. A CpG island in the promoter of the ANK1 host gene of mir-486 was shown to be highly methylated in osteosarcoma cell lines as determined by methylation-specific PCR and direct bisulfite sequencing. High methylation levels were seen for osteosarcoma patient samples, xenografts and cell lines based on quantitative methylation-specific PCR. 5-Aza-2'-deoxycytidine treatment of osteosarcoma cell lines caused induction of miR-486-5p and ANK1, indicating common epigenetic regulation in osteosarcoma cell lines. When overexpressed, miR-486-5p affected cell morphology. CONCLUSIONS: miR-486-5p represents a highly cancer relevant, epigenetically regulated miRNA in osteosarcoma, and this knowledge contributes to the understanding of osteosarcoma biology.

Rapid genome editing by CRISPR-Cas9-POLD3 fusion.

Precision CRISPR gene editing relies on the cellular homology-directed DNA repair (HDR) to introduce custom DNA sequences to target sites. The HDR editing efficiency varies between cell types and genomic sites, and the sources of this variation are incompletely understood. Here, we have studied the effect of 450 DNA repair protein-Cas9 fusions on CRISPR genome editing outcomes. We find the majority of fusions to improve precision genome editing only modestly in a locus- and cell-type specific manner. We identify Cas9-POLD3 fusion that enhances editing by speeding up the initiation of DNA repair. We conclude that while DNA repair protein fusions to Cas9 can improve HDR CRISPR editing, most need to be optimized to the cell type and genomic site, highlighting the diversity of factors contributing to locus-specific genome editing outcomes.

Comparing the epigenetic landscape in myonuclei purified with a PCM1 antibody from a fast/glycolytic and a slow/oxidative muscle.

Muscle cells have different phenotypes adapted to different usage, and can be grossly divided into fast/glycolytic and slow/oxidative types. While most muscles contain a mixture of such fiber types, we aimed at providing a genome-wide analysis of the epigenetic landscape by ChIP-Seq in two muscle extremes, the fast/glycolytic extensor digitorum longus (EDL) and slow/oxidative soleus muscles. Muscle is a heterogeneous tissue where up to 60% of the nuclei can be of a different origin. Since cellular homogeneity is critical in epigenome-wide association studies we developed a new method for purifying skeletal muscle nuclei from whole tissue, based on the nuclear envelope protein Pericentriolar material 1 (PCM1) being a specific marker for myonuclei. Using antibody labelling and a magnetic-assisted sorting approach, we were able to sort out myonuclei with 95% purity in muscles from mice, rats and humans. The sorting eliminated influence from the other cell types in the tissue and improved the myo-specific signal. A genome-wide comparison of the epigenetic landscape in EDL and soleus reflected the differences in the functional properties of the two muscles, and revealed distinct regulatory programs involving distal enhancers, including a glycolytic super-enhancer in the EDL. The two muscles were also regulated by different sets of transcription factors; e.g. in soleus, binding sites for MEF2C, NFATC2 and PPARA were enriched, while in EDL MYOD1 and SIX1 binding sites were found to be overrepresented. In addition, more novel transcription factors for muscle regulation such as members of the MAF family, ZFX and ZBTB14 were identified.

The expressed mutational landscape of microsatellite stable colorectal cancers.

BACKGROUND: Colorectal cancer is the 2nd leading cause of cancer-related deaths with few patients benefiting from biomarker-guided therapy. Mutation expression is essential for accurate interpretation of mutations as biomarkers, but surprisingly, little has been done to analyze somatic cancer mutations on the expression level. We report a large-scale analysis of allele-specific mutation expression. METHODS: Whole-exome and total RNA sequencing was performed on 137 samples from 121 microsatellite stable colorectal cancers, including multiregional samples of primary and metastatic tumors from 4 patients. Data were integrated with allele-specific resolution. Results were validated in an independent set of 241 colon cancers. Therapeutic associations were explored by pharmacogenomic profiling of 15 cell lines or patient-derived organoids. RESULTS: The median proportion of expressed mutations per tumor was 34%. Cancer-critical mutations had the highest expression frequency (gene-wise mean of 58%), independent of frequent allelic imbalance. Systematic deviation from the general pattern of expression levels according to allelic frequencies was detected, including preferential expression of mutated alleles dependent on the mutation type and target gene. Translational relevance was suggested by correlations of KRAS/NRAS or TP53 mutation expression levels with downstream oncogenic signatures (p < 0.03), overall survival among patients with stage II and III cancer (KRAS/NRAS: hazard ratio 6.1, p = 0.0070), and targeted drug sensitivity. The latter was demonstrated for EGFR and MDM2 inhibition in pre-clinical models. CONCLUSIONS: Only a subset of mutations in microsatellite stable colorectal cancers were expressed, and the "expressed mutation dose" may provide an opportunity for more fine-tuned biomarker interpretations.

Clinical and molecular implications of NAB2-STAT6 fusion variants in solitary fibrous tumour.

Solitary fibrous tumour (SFT) is a mesenchymal neoplasm characterised by pathognomonic NAB2-STAT6 gene fusions. The clinical implications and prognostic value of different fusion variants has not been clarified. In the current study, we explore the clinicopathological, prognostic and molecular differences between tumours with different fusions. Thirty-nine patients with localised, extrameningeal SFT were included, of whom 20 developed distant recurrence and 19 were without recurrence after long term follow-up. Capture-based RNA sequencing identified 12 breakpoint variants, which were categorised into two groups based on the STAT6 domain composition in the predicted chimeric proteins. Twenty-one of 34 (62%) sequenced tumours had fusions with most of the STAT6 domains intact and were classified as STAT6-Full. Thirteen tumours (38%) contained only the transactivation domain of STAT6 and were classified as STAT6-TAD. Tumours with STAT6-TAD fusions had a higher mitotic count (p=0.016) and were associated with inferior recurrence-free interval (p=0.004) and overall survival (p=0.012). Estimated 10-year recurrence-free survival was 25% for patients with STAT6-TAD tumours compared to 78% for the STAT6-Full group. Distinct transcriptional signatures between the fusion groups were identified, including higher expression of FGF2 in the STAT6-TAD group and IGF2, EGR2, PDGFRB, STAT6 and several extracellular matrix genes in STAT6-Full tumours. In summary, we demonstrate that NAB2-STAT6 fusion variants are associated with distinct clinicopathological and molecular characteristics and have prognostic significance in extrameningeal SFT.

Discovery of novel candidates for anti-liposarcoma therapies by medium-scale high-throughput drug screening.

Sarcomas are a heterogeneous group of mesenchymal orphan cancers and new treatment alternatives beyond traditional chemotherapeutic regimes are much needed. So far, tumor mutation analysis has not led to significant treatment advances, and we have attempted to bypass this limitation by performing direct drug testing of a library of 353 anti-cancer compounds that are either FDA-approved, in clinical trial, or in advanced stages of preclinical development on a panel of 13 liposarcoma cell lines. We identified and validated six drugs, targeting different mechanisms and with good efficiency across the cell lines: MLN2238 -a proteasome inhibitor, GSK2126458 -a PI3K/mTOR inhibitor, JNJ-26481585 -a histone deacetylase inhibitor, triptolide-a multi-target drug, YM155 -a survivin inhibitor, and APO866 (FK866)-a nicotinamide phosphoribosyl transferase inhibitor. GR50s for those drugs were mostly in the nanomolar range, and in many cases below 10 nM. These drugs had long-lasting effect upon drug withdrawal, limited toxicity to normal cells and good efficacy also against tumor explants. Finally, we identified potential genomic biomarkers of their efficacy. Being approved or in clinical trials, these drugs are promising candidates for liposarcoma treatment.

Integrating Anatomical, Molecular and Clinical Risk Factors in Gastrointestinal Stromal Tumor of the Stomach.

BACKGROUND: Adjuvant imatinib for 3 years is recommended to patients with high-risk gastrointestinal stromal tumor (GIST). Risk stratification is inaccurate, and risk assessments are further complicated by the increased use of neoadjuvant treatment. Anatomical criteria for prognostication have not been investigated. METHODS: Clinical, molecular, and anatomical variables were retrospectively studied in a population-based cohort of 295 patients with gastric GIST resected between 2000 and 2018. Gastric subsite was divided into the upper, middle, and lower thirds. Growth pattern was classified as luminal, exophytic, or transmural based on imaging and surgical reports. RESULTS: Of 113 tumors in the upper third of the stomach, 103 (91.2%) were KIT mutated, 7 (6.2%) were PDGFRA mutated, and 104 (92.0%) harbored genotypes sensitive to imatinib. Transmural tumors were strongly associated with a high mitotic index. Five-year recurrence-free survival (RFS) was 71% for patients with transmural tumors versus 96% with luminal or exophytic tumors (hazard ratio [HR] 8.45, 95% confidence interval [CI] 3.69-19.36; p < 0.001), and, in high-risk patients, 5-year RFS was 46% for patients with transmural tumors versus 83% with luminal or exophytic tumors (HR 4.47, 95% CI 1.71-11.66; p = 0.001). Among 134 patients with tumors > 5 cm, there were 29 recurrences. Only five patients with exophytic or luminal tumors had recurrent disease, of whom four had tumor rupture. Five-year RFS for patients with exophytic/luminal tumors >5 cm without rupture was 98%. CONCLUSIONS: In the upper third, over 90% of tumors were sensitive to imatinib. Patients with exophytic or luminal tumors without rupture, irrespective of size, had an excellent prognosis and may not benefit from adjuvant therapy.

Pembrolizumab in advanced osteosarcoma: results of a single-arm, open-label, phase 2 trial.

AIM: To evaluate the activity and safety of the PD-1 antibody pembrolizumab in adult patients with advanced osteosarcoma. MATERIAL AND METHODS: The study was a single-arm, open-label, phase 2 trial in patients with unresectable, relapsed osteosarcoma. The primary endpoint was clinical benefit rate (CBR) at 18 weeks of treatment, defined as complete response, partial response, or stable disease using RECIST v1.1. The trial had a Simon´s two-stage design, and ≥ 3 of 12 patients with clinical benefit in stage 1 were required to proceed to stage 2. The trial is registered with ClinicalTrials.gov, number NCT03013127. NanoString analysis was performed to explore tumor gene expression signatures and pathways. RESULTS: Twelve patients were enrolled and received study treatment. No patients had clinical benefit at 18 weeks of treatment, and patient enrollment was stopped after completion of stage 1. Estimated median progression-free survival was 1.7 months (95% CI 1.2-2.2). At time of data cut-off, 11 patients were deceased due to osteosarcoma. Median overall survival was 6.6 months (95% CI 3.8-9.3). No treatment-related deaths or drug-related grade 3 or 4 adverse events were observed. PD-L1 expression was positive in one of 11 evaluable tumor samples, and the positive sample was from a patient with a mixed treatment response. CONCLUSION: In this phase 2 study in advanced osteosarcoma, pembrolizumab was well-tolerated but did not show clinically significant antitumor activity. Future trials with immunomodulatory agents in osteosarcoma should explore combination strategies in patients selected based on molecular profiles associated with response.